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Coherent single electron spin oscillation in a double quantum dot system driven by a magnetic 
electron spin resonance field is studied theoretically using a Bloch-type rate equation approach. 
The oscillation frequency and relaxation time obtained using typical model parameters are consistent 
with experiment findings. The dominant decoherence mechanism is identified to be a leakage current 
through a Coulomb blockade barrier at a quantum dot during the spin manipulation. Nuclear field 
fiuctuations which induce a long relaxation time are found to contribute only negligibly to the 
decoherence despite an earlier suggestion. 
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Because of a long decoherence time, the use of elec- 
tron spins in solid-state systems, such as semiconductor 
nanostructures, is a promising approach for realizing a 
qubit, which is the basic unit in quantum computing and 
quantum information processing P, |3, Q ■ The manipu- 
lation and readout of a single electron spin arc all essen- 
tial. Solid-state implementations are expected to be most 
scalable and arc under intensive research. Electrons real- 
izing solid-state spin-based qubits are in general confined 
to semiconductor quantum dots where their numbers can 
be precisely controlled [3| . The weak coupling of electron 
spins to the rest of the solid allows long-lived spin states 
with relaxation time as long as a few milliseconds 0, Q . 
Nevertheless, this weak coupling also makes the direct 
measurement of a single electron spin challenging. Basic 
operations for building up universal spin-based quantum 
logics have only been realized recently. For example, the 
coherent manipulation of two electron spins in a double 
quantum dot (DQD) system has been achieved experi- 
mentally Q. More importantly, the coherent rotation 
of a single electron spin in a DQD system using electron 
spin resonance (_ESR) was demonstrated very recently by 
Koppens et al. [8| . Clear evidence of coherent oscillation 
of the spin state was observed. There was also a gradual 
decay of the oscillation amplitude with a time constant of 
the order of 1 /is, which was attributed to the fluctuations 
of nuclear fields due to the host materials. This is consis- 
tent with the conventional belief that nuclear fields may 
lead to fast relaxation of electron spins @, However, 
recent experiments have shown that nuclear field fluctu- 
ations can only induce the relevant spin transitions at a 
much longer time scale of 100 /is or more at large exter- 
nal magnetic fields and hence cannot be responsible 
for the experimentally observed spin relaxation. 

In this Letter, we derive a detailed theory for the ESR 
induced single spin oscillation and the associated spin- 
dependent quantum transport in a DQD system in order 
to explain the experiments in Q. By deriving and nu- 
merically integrating a set of Bloch-type rate equations 
for the reduced density matrix elements for the DQD, we 



successfully obtain the coherent oscillations of a single 
electron spin when driven by an ESR field in close agree- 
ment with experiments. We further show that nuclear 
field fluctuations cannot account for the experimentally 
observed decay of the spin oscillations in contrast to the 
suggestion by Koppens et al. 0]. It is shown to be caused 
by a leakage current through a Coulomb barrier during 
spin manipulation. Approaches for enhancing the coher- 
ence will be discussed. 

Model — The system consists of DQD connected to 
two electron reservoirs via tunneling barriers (see Fig. 1). 
We first deflne the Hamiltonian for only the DQD under 
coherent manipulation 
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where i — L ot R denotes the left or right quantum dot. 



while aL 



and 



are the creation, annihilation, 



and number operators, respectively, for electrons at dot i 
with spin a. The first four terms on the r.h.s. represent 
the electron orbital energy, interdot tunneling, and intra- 
and interdot Coulomb interactions, respectively. The last 
term describes the interactions with magnetic fields: 
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where Si is the spin operator at dot i, while g and 
are the electron g-factor and the Bohr magneton, respec- 
tively. The nuclear magnetic field at dot i due to the host 
materials is denoted by BNi. It is known to fluctuate at 
a time scale of the order of 1 s, which is much longer than 
that of the relevant electron transport process. They are 
therefore taken as stationary random fields [l^ . An ex- 
ternal field is applied in the perpendicular direction 
to generate a Zeeman splitting gfXBB^^^. Most interest- 
ingly, a spin at dot i is manipulated by the applied oscil- 
lating magnetic ESR field Bf cos{uJct) when in resonance 
with the Zeeman splitting. 
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The full Hamiltonian for the complete device of DQD 
connected to external leads is H= H^qd +-ffioads + Ht- 
The Hamiltonian of the leads is defined as i^ioads = 
Y^aka ^akaC^kcr^aka, whcrc c+j.^ {caka) IS the Creation 
(annihilation) operator of an electron with momentum 
k and spin a in lead a {a = l,r). Finally, the tunnel- 
ing coupling between the DQD and the leads is given by 

Ht = Z]feo-(^i«L(7C;fco- + ^ra^^Crka + H.c.)- 

At least one electron is always kept in the right dot 
in the experiments Q by applying appropriate gate volt- 
ages. A second electron is transported from lead I to 
r via the dots. The relevant electronic states for the 
DQD span a seven-dimensional Hilbert space. The ba- 
sis set consists of the single-electron states | 1r), \ |_r), 
and five double-electron states, namely, the double-dot 
triplets |1) = |r+) = I UU), |2) = |r_) = I UIr), and 
|3) = \To) = ^{\UIb) + \IlU)): the double-dot singlet 
|4) = \Sd) = ^(1 UIb) - I IlU)), and the single-dot 
singlet |5) = \Ss) - ^(| UIr) - I IrU))- Single-dot 
triplet states are excluded due to their much higher or- 
bital energy [ill . In this representation, i?DQD is 
rewritten as [l2| 
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+^2 coa{iJct) [ - |4) (1| + |4) (2| + H.c] , (2) 

where = ^ (Bna - Bnl ) , B, = i (Bnl + Bna) 4- B^^^t z, 
and ill. 2 = i^j^gy'siBf^ ± ^b)- have also introduced 
energy levels given by E^^ = Er ^ yf^BcrR{B^i^ + 
-^oxt)' Ei^2 = E3 ^ gfisBg, E3^4 = El + Er + Ulr, £^5 = 
2Er + Ur. a critical step in the transport is the hopping 
to or tunneling through the right dot. There can be a 
non-zero Coulomb energy barrier 



A^E5-Ei = UR- Ulr - [El - Er) 



(3) 



for the second electron at the right dot if the intradot 
repulsion Ur dominates. 

Coherent manipulation — To apply and detect the 
coherent rotation of an individual electron spin, the ex- 
periment cycles repeatedly through a manipulation stage 
and a combined readout and initialization stage (both 
stages last for 1 ^s) Q. In the combined readout and 
initialization stage, the right dot potential is pulsed low 
so that the right-dot barrier vanishes (A ~ 0). Also, the 
ESR field is turned oflF (fli = ^2 = Bf = 0). This ini- 
tializes the system at a spin blockade regime as will now 
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FIG. 1: (Color online) Schematic diagram showing electron 
transfer in a DQD structure starting at double-dot states (a) 
\UIr), (b) \UU), (c) \ UU), and (d) \UIr}. Electrons 
with 5z = in (a) and (b) admit transport under the effects of 
both nuclear fields and interdot tunneling. In contrast, trans- 
port for states with Sz = ±1 in (c) and (d) is forbidden due 
to the high PauU exclusion energy cost of the corresponding 
single-dot states. 



be explained. When a electron enters the left dot from 
the lead carrying a random spin, the DQD system takes 
one of the four double-electron double-dot states |1), |2), 
|3), or |4) with equal probability. If the system takes |4), 
the electron readily carries out a spin-independent hop 
to the right dot as A ~ and turns the electrons into 
the single-dot state |5). This transition is accoimted for 
by the Vq term in i?DQD in Eq. ([2]). The electron then 
exits to the right lead and this completes the transport 
of one electron. Alternatively, starting at state |3), the 
z components of the random nuclear fields induce a ran- 
dom relative phase between the spins. The system thus 
evolves into the near-degenerate state |4) as described by 
the B^ term in i?DQD- The electron can then be simi- 
larly transported. As a result, the second electron can 
always tunnel through the right dot for initial configu- 
rations with Sz = [Figs. l(a)-(b)]. In contrast, the 
other possible initial states |1) and |2) are split off in 
energy due to the application of a large external field 
^ext ^ V < -^N ^ and thus arc not coupled to |3) or 
1 4). They are hence spin blocked states and stop further 
transport [Figs. l(c)-(d)]. Therefore, starting from an 
arbitrary state, current flows in general until the DQD 
arrives stochastically at |1) or |2) and this completes the 
initialization stage. 

During the whole manipulation stage which follows, 
the right dot potential is pulsed up to provide a Coulomb 
blockade with A S> 0. A burst of oscillating ESR field 
is applied for a period t just before the end of the ma- 
nipulation stage (i.e., Ui,il.2 > 0). A spin which is at 
resonance can be coherently rotated. The spin blocked 
states |1) or |2) hence coherently evolves into the non- 
spin-blocked state |3) or |4) via the fli or ^2 term in 



Hbqd in Eq. ([2]). However, due to the Coulomb barrier 
A at the right dot, all transport is suggested to be com- 
pletely suppressed . We will explain later that there is 
a non-negiigiblc chance that electrons can tunnel through 
the Coulomb blockade constituting an unmeasured leak- 
age current. Other electrons from the left lead then enter 
and fill the DQD again. In all cases, there is a probability 
P(t) = P33(t) -|- P44{t) that the DQD ends up at a non- 
spin-blocked state at the end of the manipulation stage. 
Here, Pmn{T) denotes a reduced density matrix element 
for the DQD after manipulation. 

At the subsequent read-out and initialization stage, 
without the Coulomb barrier (A ~ 0), an electron can be 
transported if it is left at a non-spin-blocked state after 
manipulation. This occurs with probability P{t). Once 
an electron is transported, a random number of electrons 
may follow until a spin-blocked state is reached again as 
explained above. It is easy to show that on average one 
electron follows. As a result, the current detected in the 
readout stage is 



2e 

-frf(T) = ^[P33(T) +P44(t)], 



(4) 



where e is the electronic charge and T = 2 fis. In 
particular, at large t when coherence cannot be main- 
tained, /933(t) = P44(t) = 1/4 and the current reduces to 
/d(oo) = e/T. 

Bloch-type rate equation approach — We now derive a 
set of Bloch-type rate equations for the reduced density 
matrix elements pmn of the DQD under the coherent ma- 
nipulation by an ESR pulse. We adopt Gurvitz et a/.'s 
approach [lj| in which the many-body Schrodingcr equa- 
tion of the system is reduced to quantum rate equations 
by integrating out continuum reservoir states at the large 
voltage bias limit. We consider the high Zeeman splitting 
limit in the DQD so that spin flips caused by hyperfine 
interactions can be neglected. The rate equations for the 
diagonal density matrix elements are obtained after some 
algebra as [l3| 



Poo = -4r;poo + TrPss , 

Pmm = F; Poo - j (to I [^^DQD, p\\m) . 
P55 = -r,.p55 - i(5|[HDQD, /5]|5), 



(5) 



for m = 1,2,3 or 4 where poo = PTi?TH + PIrIr- Here, 
Ta = 2Tr Pa fl^ is the transition rate for electron tunnel- 
ing, while Pa and Qa denote, respectively, the density 
of states and transition amplitude at lead a. The rate 
equations for the off-diagonal elements are 

Cmn -p 

Pmn ^ i rPmn 

+ (1 - Cnm)r;Poo - i{m\[HuQB, p]\n), (6) 

for m,n = 1,2,3,4 or 5 (m ^ n). The coefficient 
equals one when m or n is 5 and equals zero otherwise. 
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FIG. 2: (Color online) (a) Plot of current /d(r) against ESR 
burst period r for a high Coulomb barrier A = IOOVq show- 
ing coherent oscillations; (b) Plot of measured frequency of 
oscillation f^Rabi against ESR field amplitude _B£. The solid 
line represents f^Rabi = gpBB^/2. 



Any time dependence of the coefficients in Eq. ([6]) has 
been eliminated by applying a rotating wave approxi- 
mation, which is well justified under the electron spin 
resonant condition considered here [3. 

We take typical model parameters appropriate for 
Koppcn et al.'s experiment. Specifically, the tunneling 
coupling constant between the dots is Vq ~ 0.25 peV. 
The transition rates at the leads arc F/ = = lOVb 

is chosen as 100 mT. 



17| . The perpendicular field B^^^ 
In the experiments, spins at only one dot is at resonance 
with the ESR field. We assume that this occurs at the 
left dot at a frequency lUc = gpsB^^^ with a field ampli- 
tude Bl varying from 0.74 mT to 1.5 mT. The field at 
the right dot is non-resonant and is neglected (Bfj. = 0). 
The magnitude of fluctuations of the nuclear fields was 
experimentally found to be about 2.2 mT We hence 
put = 2.2 mT and B^^^ = mT but other values of 
comparable magnitude give similar results. 

We first consider the case of complete Coulomb block- 
ade in the manipulation stage as assumed in Ref. Q 
by taking a large right-dot barrier A = IOOVq. Start- 
ing from a spin-blocked initial state, say |1), we simu- 
late the quantum transport and manipulation dynam- 
ics by numerically integrating the rate equations ([5]) and 
The current /(j(t) fiowing through the DQD as a 
function of the burst period r is then calculated using 
Eq. and plotted in Fig. 2(a). It is observed that 
/rf(r) oscillates periodically w.r.t. r. It can be well ap- 
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proximatcd by the analytical result at complete Coulomb 
blockade given by /^(r) = -^[1 - cos(riR,abiT)], 

where f^Rabi = gf^sBf^f^- This is due to a Rabi oscil- 
lation of the spin state of the electron at resonance with 
the ESR field. It hence implies an oscillation between the 
spin-blocked and unblocked states when projected onto 
the singlet-triplet basis and hence an oscillation of /^(t) 
according to Eq. ([4]). The oscillation frequency r^Rabi 
measured from our numerical data agrees with the ex- 
act relation above [Fig. 2(b)] and in particular has a lin- 
ear dependence on the amplitude of the ESR field Bf in 
agreement with experiments. In this case, we observe 
practically no decay of the current oscillations as the 
Rabi oscillations have almost constant amplitudes. This 
is however in sharp contrast to the gradual decay from 
experiments. Therefore, our result shows clearly that de- 
cay cannot be explained by the quasi-static nuclear fields 
in the presence of complete Coulomb blockade as have 
been suggested by Koppens et al. [Sll. This is in fact con- 
sistent with a recent experiment [ll| showing that the 
relevant transition between triplet-singlet spin states in 
DQD as induced by nuclear fields takes much longer du- 
rations of about 100 fis and 1 ms at external fields of 
30 mT and 150 mT, respectively. Furthermore, other 
spin-relaxation mechanisms including spin-orbit interac- 
tions and hypcrfinc interactions have even longer time 



Br = 1.5mT 



scales [201, |2l| and similarly cannot account for the decay 



of the current oscillations. 

Our theory agrees much better with experiments when 
a tunneling leakage current through the Coulomb block- 
ade barrier is considered. Here, we assume a more real- 
istic value for the Coulomb barrier A = 20Vb = 5 /xeV 
which will be further explained below. The correspond- 
ing simulation result is shown in Fig. 3. We find that the 
Rabi oscillations damp gradually during the ESR burst 
at a rate in good agreement with experiments. A smaller 
A would lead to an even faster decay. The oscillation 
period f^Rabi depends linearly on the amplitude of the 
ESR field as before. The current /^(t) oscillates around 
the limiting value /^(oo) = e/T ~ 80 fA in agreement 
with experiments. Finally, the Coulomb barrier A for an 
experimental device is limited by the fact that the right 
dot potential cannot be pulsed arbitrarily high because 
one electron has to remain all the time. From Eq. ([3]), 
one can see that A does not only depend on the differ- 
ence between the intra- and intcrdot Coulomb repulsions, 
but also on the difference of the orbital levels for the two 
dots. It can be shown that the value of A used here is 
then a reasonable one for the experiment when these two 
differences are comparable. 

In conclusion, we have investigated the coherent rota- 
tion of a single electron spin in a DQD system driven 
by an ESR field. A set of Bloch-type rate equations have 
been derived analytically and integrated numerically. We 
then obtain a current oscillation as a direct consequence 
of a coherent Rabi oscillation of an individual electron 
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FIG. 3: (Color online) Same as in Fig. 2(a), but for a small 
Coulomb barrier A = 20Vb showing coherent oscillation with 
damping. 



spin state in agreement with experiments. We also ver- 
ify that the frequency of the oscillation has a linear de- 
pendence on the amplitude of the ESR field. A leakage 
current through a Coulomb blockade barrier during the 
manipulation stage is identified as the dominant mech- 
anism of decoherence. Nuclear field fluctuations which 
induce a relaxation time much longer than the manipu- 
lation period are shown to contribute negligibly to the 
decoherence despite an earlier suggestion. Our detailed 
quantitative theory which allows the identification of the 
relevant decoherence mechanism is of great importance. 
Based on our result, it is evident that the coherence can 
be improved by suppressing the leakage current. This 
can be achieved by enhancing the Coulomb barrier A by 
fabricating a smaller right dot in a similar device, or by 
decreasing the orbital level difference via tuning the ap- 
propriate gate voltages. We hope that our results can 
motivate further experimental investigations. 
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